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Abstract 

The quantum inverse scattering method is apphed to the solution of equations for 
wave functions of compound states of n reggeized gluons in the multicolour QCD 
in a generalized leading logarithmic approximation. 
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In the leading logarithmic approximation (LLA), where g^inS ~ 1 (g' is 
QCD coupling constant), the asymptotics of amplitudes for the colourless par- 
ticle scattering at high energies is governed by the Feynman diagrams with 
two reggeized gluons in the crossing channel The intercept of the Pomeron 
trajectory turns out to be bigger than unity in disagreement with the Froissart 
bound t^l, which is a consequence of the violation of .S-channel unitarity con- 
straints in LLA. One of the possible method of unitarizing the LLA results is 
based on the effective field theory valid for the multi-Regge kinematics of pro- 
duced gluons Another (approximate) method corresponds to a generalized 
LLA, in which one calculates the contributions of the Feynman diagrams with an 
arbitrary but conserved number n of reggeized gluons in the crossing channel . In 
our previous work it was shown that the corresponding Bethe-Salpeter equations 
for the n gluon t-channel partial waves fi_j,uj = j — 1, are significantly simplified 
in the multi-colour limit N ^ oo In the impact parameter representation 
the wave functions fw{.Pii P2----Pn'i Po) for compound states of the reggeized gluons 
with two-dimensional transverse coordinates (pi, p2---Pn) have the property of the 
holomorphic factorization: 



fUPl^P2---Pn; Po) = XI /''(Pl' P2, ■■■Pn, Po)/''(Pl, P2^ -fn: Po)> (1) 

r 

where f^{f^) are holomorphic (antiholomorphic) functions of the complex coordi- 
nates Pi, p* of gluons in the impact parameter space. The position oj of singularities 
of t-channel partial waves is expressed in terms of energies e, e for corresponding 
holomorphic and antiholomorphic subsystems: 

where the total holomorphic Hamiltonian is 

n 

H = ^Hi^i+i. (3) 

i=l 

The pair Hamiltonians Hi ^ can be written in two equivalent forms 
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Hik = PrHn{pik)Pi + P^Hn{pik)Pk + ^n{PiPk) - 2V'(1) = 

(4) 

= 2in{pik) + pidn{PiPk)p-^ - 2V'(1), 

where pik = pi — Pk,Pk = = is the Euler constant. From rep- 

resentations (4) one can verify that the operator conjugated to 7Y, can be 
constructed from H, with the use of two different similarity transformations: 

i?^ = P-\..P{^HP^...P^ = p^2P23:.PnlHpiip^i...p-l (5) 

Therefore there is the following differential operator, commuting with Hamiltonian 
(3): 

A = Pl2P23-PnlPlP2-Pn, [A^H] = (6) 

and, to find a solution of Schrodinger's eqs. (2), it is enough to find all solutions 
of the differential equation : 

Af = A/. (7) 

In this paper we shall show that there is a set of n-th order mutually commuting 
differential operators, including A (6) and being a solutions of the Yang-Baxter 
equations for a monodromy matrix. 

To begin with, let us introduce the following 2x2 monodromy matrix 



^ ' » -p{d\, Q-pidij \ -pid2, 0-p2d2 



Q + Pndn, dn 
-pldn, © - pndn 



(8) 



where is the spectral parameter, dk = 
One can easily verify that its factors 



satisfy the following commutation relations 
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where the numerical matrix L,J_^f^^.j_^j^{u) is the lattice L-operator for the XXX 
Baxter model with the following nonvanishing matrix elements: 

^nai(^) = ^22^22(«) =U+1, 

^12a2(^) = ^2ljl = (11) 

1-1/2 ( . _ rl/2 _ , 
-^12;21 W — -^21;12 " 

The matrix L^/^ satisfies the three-linear Yang-Baxter equations . As a conse- 
quence of eq. (10) we obtain the commutation relations for the matrix elements 

of r(e) (8): 



T{v)i2k2T{u)i,k^Ll{l^^.^.^{u -v) = i^-/f,,fc,fc,(« - v)T{u)k2hT{v)k^n. (12) 
From these equations one can conclude in particular, that the operators 

t(G) = r(G)n + T(G)22 (13) 

commute with each others 

t{v)t{u) = t{u)t{v). (14) 
Therefore, if we expand t{Q) in the series over 

m = Qn + eQn-l+... + e'', (15) 

then the coefficients Qk, being the k-th order differential operators: 

Qk = ^ Piii2Pi2i3--Pikii^ii(^i2--(^ik, (16) 



{il,i2,...ifc} 
»l<»2<--<*fe 



commute 
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[Qfe, Qk'] — 

The operator A (7) is proportional to t(0): 



(17) 



A = i'^Qn (18) 
Another operator, which also commutes with Ti, (3) is 

Q2 = -Y.pldidk. (19) 

i<k 

Q2 is the Kazimir operator - Mi)"^ of the Mobious group of transformations: 
p concerning which eqs. (2) are invariant I^'^L 

It is plausible that all operators Qk commute with Ti. Moreover it is natural to 
expect that is a function of operators and therefore to find its eigenvalues 
it is enough to calculate eigenvalues oi Qk- Eigenvalues A and eigenfunctions of 
eq. (7) can be obtained by using the algebraic Bethe ansatz t^l. In the future 
publications we shall use the quantum inverse scattering method, developed by L. 
Faddeev and his collaborators to solve eqs. (2), (7), starting from commutation 
relations (12). 

The remarkable mathematical properties of the equations for QCD scattering am- 
plitudes at high energies could be explained from the fact, that the Yang-Mills 
theory is a low energy limit of the super-string model, having a very rich symmetry 
group. 

The author is thankful to the Dipartimento di Fisica deU'Universita di Padova for 
the hospitality and to Prof. A. Bassetto and Dr. M. Matone for helpful discussions. 
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